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Abstract 

An extension of a model of neutrino masses to the quark sector provides an 

interesting link between these two sectors. A parameter which is important to 

describe neutrino oscillations and masses is found to be a crucial one appearing 

in various "penguin" operators, in particular the so-called Z penguin. This 

parameter is severely constrained by the rare decay process Kl ^^/^^ • This 

in turn has interesting implications on the decay rates of other rare processes 

such as Kl — > ne, etc., as well as on the masses of the neutrinos and the 

masses of the vector-like quarks and leptons which appear in our model. 
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I. INTRODUCTION 



In the last few years, we have witnessed a flurry of far-reaching experimental results, 
among which are the neutrino oscillation data [|^, data on direct CP- violation such as e'/e 
and upper bounds on Flavour- Changing-Neutral- Current (FCNC) rare decays of the 
kaons. On the one hand, the neutrino oscillation data clearly points to possible physics 
beyond the Standard Model (SM). On the other hand, it is still not clear if the new results 
on e'/e, which differ roughly by a factor of two from present calculations within the SM, 
imply any new physics since the aforementioned calculations are still plagued with non- 
perturbative uncertainties. For the kaon's FCNC rare decays, the experimental situation 
is still far from giving evidences of physics beyond the SM or to confirm the SM itself. 
Nevertheless, whatever new physics, which might be responsible for giving rise to neutrino 
masses, could, in principle, affect the quark sector, and hence, also quantities such as e'/e 
or the branching ratios of the kaon's FCNC rare decays . If this is the case, results from 
the quark sector could then be used to put constraints on the lepton sector itself since it 
is possible that both sectors have a common set of parameters, a desirable feature of any 
model which purports to deal with issues of fermion masses. The aim of this paper is to 
explore this possible connection between the two sectors. 

The subject of neutrino masses has been invigorated in the last few years due to new 
results on neutrino oscillations which suggested the possibility of a non- vanishing mass for the 
neutrinos. Models have been (and are still being) built to try to describe these oscillations. 
In a large number of cases, efforts were mainly concentrated on the type of neutrino mass 
matrices which could "explain" the oscillation data, with very little attempt made in trying 
to connect that kind of physics to the hadronic sector. However, it is perfectly reasonable 
to expect that the two sectors are somehow deeply connected, and that a constraint from 
one sector can give rise to constraints on the other sector. This particular connection will 
be the subject of the present manuscript. This paper is an extension to the quark sector of 
a model of neutrino mass presented in Ref. [Q. In consequence, we shall show that there is 
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a set of parameters which appear in both sectors, and that the constraints obtained in the 
quark sector have interesting imphcations on the neutrino sector. In particular, we shall use 
recent results from limits on various rare decays to constrain a common parameter which 
appears in both sectors. ( As we will show, because this common parameter is real, there 
will be no new contributions from our model to e'/e.) 

The plan of this paper is as follows. First, we present our model for the quark sector 
based on a previous model for the neutrino sector |^. We then proceed to enumerate and 
compute various FCNC operators which arise in this model. These operators are important 
in the analysis of various rare decays. Finally, we will use these FCNC operators in the 
computation of their contributions to the aforementioned quantities, and set constraints on 
the parameters of the model. Using these constraints, we look at the question of how they 
would affect the neutrino sector. We will also show at the end of the paper that there are 
some interesting correlations between the value of the branching ratio for Kl ^ fie and 
the mass of the weak-singlet quarks and charged leptons which appear in our model. When 
these particles are "light" enough to be produced at future accelerators, the branching ratio 
for Kl ^ fie is too small, of order 10~^^, to be detected, while for a branching ratio which 
could conceivably be measured in a not-too-distant future, e.g. of order 10~^^, the masses 
of these singlet quarks and leptons will be too large, a few hundreds of TeVs or more, to be 
produced by earthbound laboratories. 

II. EFFECTIVE VERTICES FOR FCNC PROCESSES IN THE MODEL OF 
REF. [§ 

The model used in Ref. [0] to describe the lepton sector is summarized in Appendix A. 
(A quick look at that appendix will help with the notations and particle content.) In this 
section, we will simply write down the part of the Lagrangian for the quark sector which 
is relevant to the construction of various FCNC operators. First of all, the common link 
between the quark and lepton sectors is in the scalar sector. These are the SM Higgs field 
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(f) and the Family symmetry Higgs field Q". As we shall see below, it is these scalars which 
"transmit information" from the quark to the lepton sector and vice versa. In particular, as 
can be seen from Appendix A, Q°' provides the common set of parameters which appears in 
the two sectors. The generahzation from the lepton sector to the quark sector necessitates 
the introduction of a new set of vector-like quarks, F*, M\, and Ml, where i is a color 
index, in perfect analogy with the leptonic vector-like fermions, F, J\4i, and M.2- As we 
shall see below, these vector-like quarks can have masses as low as a couple of hundred 
GeV's, making them very attractive for potential discoveries at the LHC. We shall come 
back to this interesting issue in the last section of the manuscript. 

The particle content and quantum numbers of the model are listed in Table 1. Notice 
that, among the scalars listed in that table, there is one which does not have the appropriate 
quantum numbers to be able to couple directly to the quarks: p". The gauge structure of 
the model (for both quarks and leptons) is given by: 

SU{3), ^ SU{2)l ® U{1)y ® SO{A) ® SU{2),^ , (1) 

where SU{2)^^ applies only to the right-handed neutrinos, and 5'(9(4) is the family symme- 
try. The form of the quark Yukawa Lagrangian is similar to the one from the lepton sector 
(Appendix A), with the introduction of the following new parameters: 

gu, gd, Gi, G2, Gs, Gmi,Gm2^ ^F, A^Mi, AfA42 • (2) 

The quark Yukawa Lagrangian takes the form 

^quark = gdql(t>daR + guOl^UaR + 

GiqI^o^Fr + Gm.FlcI^Mir + GM,FLm2R + 
MfFlFr + MmMilMir + 

MmM2lM2r + h.c. , (3) 

with the important difference with respect to the lepton sector in that there is no coupling 
between the quarks and the scalar before symmetry breaking. After symmetry breaking, 
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the mass eigenstates are linear combinations of Q and p, and vice versa, as shown below and 
in Appendix A. 

The way the FCNC processes take place in the quark sector is identical to the way they 
operate in the lepton sector i.e. via loop diagrams. In the SM, as it is well known, the flavour 
diagonal structure of the basic vertices involving 7, Z and g forbids the appearance of FCNC 
processes at the tree level. However the FCNC processes can happen at one loop or higher 
order level, mediated by a combination of flavour changing charged currents coupled to the 
W's. The fact that these processes take place only as loop effects makes them particularly 
useful for testing the quantum structure of the theory and in search for physics beyond the 
SM. In our model, beside the ways in which FCNC processes can happen in the SM, it can 
also happen because of the couplings of quarks of different flavour to the same vector-like 
fermion F and to the NG bosons f2j and the pseudo NG bosons Repi. This is made possible 
by the mixings among the NG bosons Qi and the pseudo NG bosons Repj with different 
family indices i and j. We denote the relevant scalar mass eigenstates by Cli and Repi , in 
terms of which we can express the states entering the Yukawa Lagrangian by 



The states fij are the NG bosons which are absorbed by the corresponding family gauge 
bosons. When the NG bosons get mixed, there will be mass mixings among the correspond- 
ing family gauge bosons. If we denote by Aq the orthogonal matrix which diagonalizes the 
family gauge bosons mass matrix (Ref. 0), we can express the states Q'^, corresponding 
to the longitudinal components of the gauge boson mass eigenstates, in terms of the mass 
eigenstates Qi 



fij = cos — sin (3Repi 



(4) 



Repi = sin + cos l3Repi . 



(5) 




(6) 



with Aq given by 
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(7) 



This mechanism, introduced to describe neutrino mass sphtting in the neutral lepton 
sector, gives the following terms in the Yukawa Lagrangian: 



dgiai^R = cos(3G,qin,Fn - sin PGiq^^RekFR 

= cos(3Giq[Al^ijQ'^Fji - sm(3GiqiAj,^,^Rep'^FR , 



(8) 



G2MiL^id)^ = cos (3G2MiL^idi^ - sin/^GsA^iL^epid'^ 

= cos (3G2MiLA^ ijQ'jd'ji - sin (3G2MiLAf^ ijRepjd'^ , 



(9) 



where we assume that the same Aq diagonalizes also the mass matrix of the pseudo NG Repi 
(Ref. [0]). It is important to notice here that H4 and h^, whose mass eigenstates are given 
by H4 and h^, see Eq. 0, are not coupled to the quarks of the first three families. This 
implies that they will not propagate in the loops of the diagrams describing processes with 
external quarks of the first three generations. In this paper we will look only at processes 
involving the quarks of the first three families and so we will not care about the presence in 
our model of H4 and /14, where 



H4 = cos a H4 + sin a , 
= — sin a H4 + cos a . 



(10) 



In the following we will present the expressions of the effective vertices for FCNC pro- 
cesses mediated by Z, 7, and g. There is no extra effective vertex with the SM Higgs in our 
model because the vector-like fermions do not couple to the SM Higgs field. We will also 
compare these expressions with the corresponding ones from the SM. 

For all the effective vertices, the expressions will be given in terms of the linear combi- 
nation 
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M^Mncos^P + MpSm'^P, (11) 

where AIq and Mp are the contributions to the effective vertices when the particles prop- 
agating in the loops are respectively the NG bosons Q\ and the pseudo NG bosons Rep[. 
From the lepton sector we have estimated tan (3 to be 

t^H^-^sl-^, (12) 

with gp being the 50(4) gauge coupling. V and V are the vacuum expectation values (VEV) 
of p and Q respectively, Mp and the masses of the vector-hke fermions introduced in the 
lepton sector, and Mq the central value for the masses of the family gauge bosons. Taking 
gp ~ 0(1), and, for the masses, the values required in the lepton sector to have a proper 
value for the neutrino of the 4*'* generation, tan /3 turns out to be much smaller than unity. 
This makes the contribution to M. due to the pseudo NG bosons neghgible, being suppresed 
by the factor sin^ (3. In the following we will give only the expressions for the contributions 
due to the NG bosons Q\. All the expressions that will be given for the effective vertices will 
have to be multiplied by the factor cos^ /3. We will focus on the expressions for the effective 
vertices in the transition s — > d. 



III. THE Z EFFECTIVE VERTEX 

To calculate the Z effective vertex we have to sum up all the contributions coming 
from the different vector-like fermions propagating in the loop diagrams. We will present 
separately the expressions for the amphtude due respectively to the Fa and M.i vector-hke 
fermions. We remark here that Fd is the down component of a doublet of SU{2)l while Mi 
is a singlet of SU{2)l ( see Table 1). Because the dominant terms of the amplitudes are 
proportional to T^l of the vector-like fermions, the contribution due to F^ is singled out. 
The contribution from M.i will be suppressed by the factor I/Mq with Mq being the scale 
of breaking of the family symmetry. The amplitude due to F^ is as follows 
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(13) 
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9+ 



-9A 



-T- 



3L 



(14) 



It is important to notice here that with Mq. we have indicated the poles in the propagators 
of the longitudinal components of the gauge boson mass eigenstates fl[. 

Using for Mq. the eigenvalues of the family gauge boson mass matrix in the simplistic 
form introduced for the lepton sector 



Ml = Ml 



b 1 
1 



(15) 



given by 



Ml = M2(l+6), 
Ml = M2(l-6), 



Ml = Ml, 



(16) 



where 6 is a small parameter less than unity, and the orthogonal matrix is given in Eq. 
1^, we obtain the following final expression for Ai^^ 



M 



F,fi _ ~'^9 „2 

^ IGvr^ cos(6'vy 



9'pCo{xF,b)-sY{l-r)d 



(17) 



with 



Co{xF,b) 



A{1 + b - xpf i-l + b + xpf 



-l + bfil + b-xpf ln(l-6) 



+ (l + 6)'(-l + 6 + x^)' ln(l + 6) 

-2bxF if - {-l + xpf + 2 (-l + b'^ + xf) ln(xF 



(18) 



where xp is defined as 



Ml 
Ml 



(19) 



Making the Taylor expansion of CQ^xpih) in the parameter h we found that the first non 
zero contribution comes from the hnear term in h. So we can rewrite C(){xf-i h) as 



Co(xF,6) =6CoM + C»(6'), 



(20) 



with Cq{xf) given by 



Cq{xf) = - 



\ + Xf 



+ 



2x, 



1-xfY {I-xfY 



In xf 



(21) 



M-z ^o"^ takes the simple form 



^ 167r2 cos(^w^) 



glbCo{xF)sril-l')d+0{b^), 



(22) 



which shows, in the first approximation, an explicit linear dependence on the parameter b. 
In Appendix B, it will be shown in a simple example the mechanism that produces the b 
dependence of the amphtude for FCNC processes and it will be compared with the GIM 
mechanism. 

The contribution to the Z effective vertex due to Mi takes the form 



where q is the momentum transfer and with 

'b (-5 b^ + (-1 + XMf (5 + xm) (-1 + 7xm)) 



(23) 



Do{xM,b) = 



{l + b-xufi-'^ + b + XMf 
2b^ {5 + XM {~22 + 23xm)) 
(1 + 6 - xm)^ (-1 + b + xm)^ 
3 (-1+6)' (-1 + & + 3xm) ln(f 

(-1 + 6 + 2;^) 
3(1 + 6)' {l + b-3xM) Hf- 



XM 



M ] 

1+b) 



[1 + b- Xm) 



Xm_ 

27 ' 



(24) 



where xm is defined as 



Ml 



(25) 



It is important to notice here that the structure of the operator in Ai'^^'^ is similar in 
character to that for the photon. In fact, = 0. This happens because M.i is 

a singlet of SU{2) and so the current of A^i that couples to Z has the vectorial nature 
characteristic of the currents interacting with the photon. Making a Taylor expansion of the 
function Dq{xm, b) in the parameter b we found also that in this case the first non vanishing 
term in the expansion is that one linear in b. So we can rewrite Do{xm, b) as 



Do{xM,b) = bDo{xM) + 0{b^), 



(26) 



with Dq{xm) given by 



Do{xm) 



Now Ai^^'^ takes the simple form 



27 



11 + 61nxM — 18x m^^xm — 63xm + 45x|^ + 7x\^ — 36x\j\nxM 

(I-xm)" , 



(27) 



^ 167r2cosf^w)4M2 



7w) ^^^Mi 



■bDoixM) s(gV - + 7')^^ + 0{b^) . 



(28) 



For comparison we show the explicit expression of the corresponding contribution from 
the SM, as given in Ref. H. 



M 



SM,ii 



gX:VuCl^\xt)-si^{l-^')d.^ 



IGvr^ cos( 



(29) 



where 



Co'^ixt) = V 



Xt 



Xt-Q ^ 3xt + 2 
xt-l [Xt - ly 



(30) 



with Xt defined as 



Xt 



(31) 



and g is the SU{2)l gauge coupling. 
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IV. THE 7 EFFECTIVE VERTEX 



The effective vertex of 7 comes also from the contributions due the and Aii vector-hke 
fermions. What is remarkable is that these two contributions will contain the same function 
Dq (Eq. |2^), which appears in the contribution of A4i to the Z penguin. This happens 
because, as mentioned before, A^i is a singlet of SU{2) and is coupled to Z in the same way 
as to 7 except for the coupling constant. The contributions to the effective 7 vertex from 
Fd and Aii will differ from each other in the form of the operators and with xp and xm 
appearing in Dq respectively. We obtain the following final expressions 

= T^^iii " ^'^^^^ ~ ^'^^ + ^^^'^ ' ^^^^ 

and 

M^-^ = ^,j§^bDQ{xM)s{q'l' - q'm + l')d + Oib') . (33) 
The corresponding expression for the SM as given in Ref. M is 



-M^" = t£~2^KVuD',''{x,) rsiq'r - - l')d. (34) 



with 



D^^ixt) = — lnxt + + * ^ * — — ^Inxf. (35) 



V. THE GLUON EFFECTIVE VERTEX 

The function Dq which appears in the 7 vertex also appears in the gluon effective vertex, 
with the only difference appearing in the prefactor due to the different nature of the charge. 
This is so because g and 7 have vector couplings to the same particles in our model. This 
does not happen in the SM where 7 can couple to W while the gluon cannot. As a result, 
the number of diagrams in the SM, used to describe the effective vertices, and the functions 
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that appear in there, are different for the gluon and for 7. We introduce the function Eq 
given by 

Eo = -3Do , (36) 

to absorb the factor —1/3 in Dq coming from the electric charge of the quarks. The final 
expressions for the contributions to the gluon effective vertex due to Fd and A^i are (a, (3 
are color indices and gs is the SU{3)c coupling) 



and 

= l^ii^^^"^"^^^'"^^'^' " ^'^^^^ + + ^^^'^ ' ^^^^ 

where (T^^)" is the a-th generator of SU{3)c. The corresponding expression for the SM as 
given in Ref. |Q is 

M2^sM = iit^^^t:vt,E^'\x,)Wi'r - q'm - i')n,dp , (39) 

with 

£;5M. ^_ 2 x,^(15-16x, + 4xg) x,(18-llx,-x,^) 

Knowing the expressions for the effective vertices in the new model is what we need to derive 
completely the new contributions to the FCNC processes with AS" = 1 involving penguin 
diagrams. These new contributions differ from the corresponding ones in the SM only in the 
effective vertex. 



VI. BOX AMPLITUDES 

For the contributions to the FCNC processes with AS = 1 coming from box diagrams 
we have to calculate directly the amplitude. In Fig. 5 we show the new one loop diagrams 
describing FCNC processes with AS* = 1 coming from our model. The new contributions 
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coming from the box diagrams to the FCNC processes with AS = 1 are suppresed by the 
factor 1/Mq. This makes these box contributions neghgible with respect to the dominant 
contribution, due to the Z penguin, which is suppressed only by the factor coming 
from the Z propagator. The number of box diagrams is eight, considering only the ones in 
which the NG bosons propagating in the loop are the The diagrams in which the pseudo 
NG bosons Rep[ are propagating in the loop are suppressed by sin^/? while those which 
contain both NG bosons and pseudo NG bosons are suppressed by sin^ p. The first two 
diagrams (Fig. 5a,b) have the F vector-like particles propagating in the loop. Neglecting 
terms of order m^/M^ with rriq being the mass of the external quarks, we find the same 
amplitude for the two diagrams. We also use the identity 

s{p^)r{l - l')d{p,)d{p,)Y{l - l')dip2) , (41) 

applying Fierz transformations. For both diagrams the amplitude is given by 

Mto. = jMrri Boi^F, b) 57^(1 - 7')^ d^,{l - ^')d , (42) 

where Bq{xf, h) is given by 

-(1 + hf + x'^ + 2 {l + h) X (ln(l + 6) - ln(x)) 



2 1 + — X 



;i -hf + x'^ + 2 {l-h) X (ln(l - 6) - ln(a;)) 
2(1 -6 -x)^ 



(43) 



For the function BQ{xF,b), in the same way as for CQ{xF,b) and DQ{xF,b), we make a 
Taylor expansion in the parameter b. Also in this case as in the case for the Z and 7 vertices 
the zeroth order term in the expansion is not present and the first non vanishing term is 
that one which is linear in b. The amplitude in Eq. ^ now takes the simple form 

M,o. = Y^^^Bo{xf) -sY{l - l')ddr{l - l')d + 0{b') , (44) 

G 

where Bq{xf) is given by 
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^, , l + {A-5x) x + 2x {2 + x)\n{x) 

M^f) = ^ , , , ,4 ■ (45) 

(-1 + X) 

The amplitudes for the box diagrams in which we substitute the vector-hke fermions 
F with Ail have similar expressions except in the operator which now has the form {V + 
A){V + A) instead of (V — A){V — A) and with xm now appearing in the function Bq. As 
in Ref. we expect M_M-^to be much larger than Mp so that xm = M^_^/Mq is closer 
to unity while xp = M'^/Mq can be much less than unity. As can be seen from Fig. 6, 
Bq varies by at most an order of magnitude for < xp^M < 1- As a result we expect the 
contributions of the two diagrams due to A^i to be smaller than the corresponding ones due 
to F by not more than one order of magnitude. For the four box diagrams in which there 
are one vector-like particle of type F and one of type A^i, the analytical expressions are too 
complicated to be written down in this paper. Instead we perform a numerical evaluation 
for the function Bmixo{xp,XM) in the mixed case giving us values between Bq{xf) and 
Bo{xm)- The operators for these last four box diagrams are of the form {V — A){V + A) 
and {V + A){V - A). 

The amplitude for the box diagrams in the SM is given by 

-Mfox = V*Vu\V,,\' B^,'\x,) .-7^(1 - ^'>)ddY(l - f)d , (46) 

where 



Xt ^ Xt In Xt 



l-Xt {xt - 1)2 



(47) 



VII. CONSTRAINT ON b FROM THE UPPER BOUND ON Kl 

From our analysis of the new contributions to FCNC processes with AS* = 1, it turns 
out that the main contribution is given by the Z penguin, all the other contributions being 
I/Mq suppressed. 

The FCNC processes with AS* = 1 in kaon physics are reasonably well described by 
the SM, with uncertainties coming from non perturbative QCD effects. On the other hand 
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the apparent discrepancy for example between the SM estimates and the data invites for 
speculations about non standard contributions to e'/e 0, and also for FCNC rare deacys 
there are still margins for effects of new physics p[. Now in our extension of the SM, the 
only non negligible contribution comes from the Z penguin and corresponds to the operator 
V — A, the same as in the Z penguin of the SM. 

The way we proceed follows Ref. [Q. First of all we write down the effective La- 
grangians corresponding to the flavour-changing coupling of the Z boson to down-type 
quarks in the SM and in our model, looking only, for the latter case, at the dominant 
contribution. We have 

= 1,2%, , /^r^'7"(l - l')dZ, + h.c. , (48) 

with 

Zi"" = V*VuCt'{x,) , (49) 



and 



= , \a Zlm^ - l')d Z, + h.c. + Oib') , (50) 

XOTT COS 



with 



Z!, = bCo{xp). (51) 



We notice here that Eq. and 50 are related to Eq. and through the equation 



C = -iM^Z^ + h.c. . (52) 

As showed in Ref. [|] the coupling Z^^ is complex, being the product of the quantity V^* V^d, 
which is complex, where Vij are elements of the CKM matrix, and the function Co(xt) which 
is real. In particular, from the standard analysis of the unitary triangle, one has 

Im{V*ytd) = (1.38 ± 0.33) ■ 10"^ , (53) 
Re{VZVta) = -(3.2 ± 0.9) ■ 10"^ , (54) 

15 



which, using for Co{xt) the value of 0.79 corresponding to the central value of the top quark 

mass, fht{mt) = 166 GeV, give 



;i.09±0.26) ■ 10-\ 
-(2.54 ±0.71) ■ 10-^ 



(55) 
(56) 



The complex nature of the coupling Z^^^ is responsible for CP violation in the SM, as can 
be verified by looking at the expression for e'/e as given in Ref. P]: 



e 
e 



(57) 



Rest 



where 





'l.2-R,\rP\Bi'/''>' 







(58) 



and 



e 



td 



Rest 



-2.3 + Rs 



l.l\rf\B'i"^ + {l.Q + Q.l2\rf\)Bf"^ 



(59) 



which is proportional to Im{y^gVtd)■ AH the parameters appearing in the above expressions 
are fully described in Ref. 0. As pointed out in Ref. if we assume that no new operators 
in addition to those present in the SM contribute, and this is true in the approximation of 
neglecting the new contributions which are suppressed by the factor l/M^, the replacement 



Zds, where 



SM I 9f yF 
2 '^ds ' 



(60) 



a relation that holds separately for ReZ^s and ImZds 



ReZ, 



ds 



ReZg'^ + %ReZ^ 
9 



ds 1 



ImZ,, = ImZ',f + %ImZ^ 
9 



ds 1 



(61) 



(62) 



which is justified without the modification of QCD renormalization group effects evaluated 
at NLO level for scales below 0{mt). This means that to look at the effects of new physics, 
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described by a modified effective coupling of the Z boson to down-type quarks, to the 
quantity e'/e, we just need to substitute ImZ^s to ImZ^^^ in Eq. Looking at Eq. |51|, 
Z^g is the product of the parameter 6, introduced in Eq. |T^, which has been chosen real 
and the function Co{xf)i which is also real making Z£ real. This implies that in our model 
there are no corrections to ImZ^^ and consequently to e'/e. Obviously this means that we 
cannot use the experimental results on e'/e to constrain the parameter h. As shown in Ref. 

if new physics affects ReZds, as it is the case in our model, the process to look at to 
reveal effects of new physics is the FCNC decay Kl — ^ ^~^fi~ , whose experimental branching 
ratio BR{Kl = (7.2 ± 0.5) ■ 10^^ |jl2|. The effects of new physics appear in the 

short- distance (SD) contribution to BR{Kl /U^/U") 

BR{Kl fi^fi-)sD = 6.32 • 10-3 [ReZds - BoRe{VZVu) + A,] ' , (63) 

where Bq = —0.182 is the box diagram function evaluated at Tht{mt) = 16QGeV, and 

Ae = -(6.54 ±0.60) ■ 10-^ (64) 



represents the charm contribution []TT|. From the analysis in Ref. |^ of the long distance (LD) 
and SD contributions to BR[Kl — > /x^/i") the highest possible value for BR{Kl — > fJ.^fi^)sD 
is derived 

BRiKL ^ fi^li-)sD < 2.8 • 10-^ . (65) 



Now using Eq. |6^ together with Eq. ^ ^ and ^ we obtain the following upper bound on 
the parameter b 

6<2.5-10-\ (66) 

where in Z£ we have used Co(10-^) ^ 0.5 which corresponds to Mp = 200Ge\^ and Mq = 
200Tel^, and in Z^s we have used g = 0.65 and gp = 1- 

Looking at the expressions for BR{Kl — > tt^uu), BR{Kl vr^e^e^) and BR{K^ — >■ 
TT+z/z/) given in Ref. f^, we observe that there are no new contributions from our model 
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to the first two decays, which depend only on ImZds, while there is a contribution to 
—>■ TT^uu, which depends also on ReZ^s- In Ref. ||^, the expression for the upper bound 
on BR{K^ —>■ n^uu) is given in terms of BR{Kl tt^uu), and the parameter k related to 
BR{Kl fJ'^fi''~)sD by the relation 

BR{Kl ^ /iV~)5D = ■ 10-' . (67) 

Now, if as in Ref. p we use the upper bound for BR{Kl fJ'^fi'^)sD given in Eq. ^ we 
obtain 

BR{K+ TT+z/p) < 0.229 ■ BR{Kl n^i^u) + 1.76 ■ 10~^° , (68) 

while, if we use the upper theoretical value for the BR{Kl — > Ai+/i^)||f given in Ref. 0, 
we obtain 

BR{K+ TT+z/p) < 0.229 ■ BR{Kl ^ 7r°z/z/) + 1.15 ■ 10"^° . (69) 

In Ref. 0, for the particular scenario in which all the effects of new physics are encoded 
in the effective coupling Z^si BR{Kl —>■ tt^uu) is estimated to range in the interval (1.3 • 
10^^*^, 2.4 ■ 10^^°) for B^^"^^ = 0.6, and this makes the contribution due to the maximum 
value for BR{Kl fJ'^fJ'~)sD in the expression for the upper bound of BR{K^ n^uu) the 
dominant contribution. The difference between the two upper bounds for BR{K^ n^uu) 
shows the room for a possible contribution of new physics associated to the quantity ReZ^s- 
Using the constraint on b derived in Eq. ^ we can now check if the amplitudes that we 
know to be 1/Mq suppressed are negligible with respect to the corresponding SM amplitudes. 
In the case of the photon, matching the absolute value of the functions which multiplies the 
V — A operator for the new amplitude and the SM one, we have, after cancelling out all the 
common factors, 



Ml 



g^pbDoixp) 



<1.7-10-\ (70) 



gme{V,lVtd)D^^\xt) 

where we have used Do{lO'^) ~ -2.66 -10^^ and for the SM quantities D^^^{4:.27) -0.46, 
Re{V^*Vtd) as given in Eq. and obviously the constraint on b from Eq. 



In the case of the gluon we have 



Ml 



g],bEo{xF) 



<8.7■10-^ (71) 



where we have used Eq^ (4.27) ^ 0.27. From what we have seen above, the amplitudes 
with the V — A operator for the photon and for the gluon, which have corresponding ones 
in the SM, are neghgible. Now the amphtudes with the V + A operator which appears in 
our model, and are present in the effective vertices for the Z, the photon and the gluon, are 
also 1/Mq suppressed and are not relevant. 

For the box diagrams, we look at the function which multiplies the {V — A){V — A) 
operator. We obtain the ratio 

4gj,bBo{xF) 



Ml 



<0.19, (72) 



g^Re{y^lVtd)\Vu\''Bl^'{xt) 

where we have used -Bo(10^^) — 1 and for the SM quantities Bq^'^{A.27) ~ —0.18 and 
\Vtd\ — 9.1 ■ 10^'^. Similarly, the amplitudes containing the operators {V + A)(y + A), 
(y — A)(y + A) and (V + A){y — A) are also not relevant. It is important to notice here 
that, for the box diagrams, the 1/M^ suppression of the new contribution is balanced by 
the \Vtd\^ suppression of the SM one. 

In Ref. it has been shown in two numerical examples how it is possible to derive the 
neutrino mass splittings Amgg and ^m^^ for different values of h. In fact, in Ref. 0, h has 
been introduced as a free parameter chosen to be smaller than unity. In this paper it has 
been shown how it is possible to constrain the parameter h by looking at the quark sector. 
In accordance with the upper bound for h derived in this paper and presented in Eq. |6^, 
we give here the neutrino mass splittings obtained numerically in Ref. for the particular 
choice of the parameter h = 0.000095 

Am^2 = 2.02 ■ IQ-^eV^ , (73) 
Aml^ = 5.497 ■ IQ-^eV^ . (74) 
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It is important to notice here that, when one chooses a particular value of the parameter 
b, one also has to choose the internal loop variables which appear in the function given in 



Eq. |92|, in order to reproduce the mass splittings of Eq. |73| and |7^. This means, as one 
can see from the two numerical examples with different values of the parameter b in Ref. 
0, that different values of b will correspond to different sets of masses for the three light 
neutrinos. 



VIII. OTHER RARE KAON DECAYS BEYOND THE SM 

Being sensitive to flavour dynamics from few MeV up to several TeV, rare kaon decays 
provide a powerful tool to test the SM and to search for new physics. Decays like Kl —>■ fie 
and Kl vr/xe are completely forbidden within the SM p], where lepton flavour is conserved, 
and are also absolutely negligible if we simply extend the model by including Dirac type 
neutrino masses with the standard Yukawa mass term. 

In our model the neutrinos are still only Dirac, but the way they get masses is through 
loop diagrams, and processes like Kl —>■ yue or Kl vr/ie are made possible through the 
exchange of virtual NG bosons They, in fact, as we have already said, can couple to 
different flavours and, in the way we build our model, are the same for the lepton and the 
quark sectors. 

In the following, we will calculate explicitly the branching ratio for the decay Kl fie 
in our model, and we will give a theoretical range of values corresponding to different values 
for the mass of the vector-like fermions Mj^i,q. 

This decay happens through five diagrams, shown in Fig. 8: four box diagrams that have 
to be considered with their corresponding crossed diagrams, and the diagram obtained by 
linking two effective Z verteces by a virtual Z. It can be shown that the crossed diagrams 
are b suppressed with respect to the ones shown in Fig. 8. The diagram where the virtual Z 
is exchanged is also b suppressed, but not 1/Mq suppressed as the surviving box diagrams 
and so we will expect the sum of all the contributions will depend on the values of the 
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parameters b and Mg- The final amplitude for the process Kl — ^ fie is given by 



M 



'9l 
47r2 



327r4M| cos2 (Ow) 



+^B{xMg , a^F, , b)s'y'' Pr d 67^, Pl /i 
^2 



(75) 



with P being the function coming from the box diagrams (which is too long to be written 
down explicitely in this paper), Cq given in Eq. ^ and Pl^ = To obtain the BR we 

can use the expression given in ref. ^ for our case 

'gp lOOTel/'^ 



-12 



BR= 11.24-10 

9 Mg 

that for our purpose can be rewritten as: 



{Clq - CRq) {C^i + C^i) 



(76) 



BR = 11.24- 10 



-12 



gp lOOTel/' 



[{CLqCuY' + {ClqCriY + {CRqCaY + {Cr^CriY 



9 Mg 

—'^{CLqCLl){CRqCLl) — '^{C LqC Rl){C RqC Ri)] , 

where C^g, CRq, Cli, Cri appear in the expression for the operator Ov,a 



(77) 



V,A 



9f - 



2Ml 



sYlCLqPL + CRgPR]de^^[CLlPL + CRiPRjfl + h.c. . 



(78) 



From the way we have written the amplitude in Eq. ^ it is possible to isolate the four 
contributions CiqCu^ CLqCm, CRqCu and CRqCRi that appear in Eq. ^ 



- ^'^^^^ -b'c,{xF^)c,{xF:, 



CLqCu = -^B{xFg,XFi,b) - 
ChqCRl = ■j^B{xFg,XMi,b) , 



327r4M| cos2 



47r2 
^1 



(79) 



CRqCu = -—^B{xMq,XFi,b) 



CrinCh 



An 
9l 



B{xM„,XMi,b) . 



^Rq^Rl ^^2^ 

In Fig. 9, the logarithm of the quantity (C^g — CRq^iCli + C'ri) is plotted as a function 
of the logarithm of the ratio xm/xf-, with Mg = 200Te\^, xp = 10~^ and using for b the 
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upper bound given in Eq. so that the contribution coming from the diagram with the 
virtual Z has its highest value. In doing the plot of Fig. 9 we have also chosen Mj^i = Mj^i^ 
and Mpi = Mp, = 200^6^. In Fig. 10 the logarithm of BR{Kl fie) is plotted as a 
function of the logarithm of the ratio xm/xf, for the same values of the quantities Mq, xp 
and b as in Fig. 9, and the same choice Mj^^i^ = Mj^q and Mpi = Mpq = 200GeV. It can be 
seen from Fig. 10 that BR{Kl —* fie) depends strongly on xm, starting from 6.65 ■ 10^^^ for 
xm/xf = 1, and reaching the asymptotic value of 1.02 ■ 10~^^. It is important to notice here 
that the lowest value for BR{Kl — > fie), when xm/xf = 1, corresponds to the situation in 
which the contributions of the box diagrams cancel out and the only contribution left is due 
to the diagram with the virtual Z. In this case we have 

BR{Kl fie) = 11.24 ■ 10" 

= 6.65■10"2^ (80) 

where we have used for b the upper bound given in Eq. ^ The asymptotic value for 
BR{Kf — > fie), when xm/xf ^ 1, corresponds instead to the situation in which the box 
diagram with and F'^ propagating in the loop, gives the dominant contribution, and 
all the contributions from the other diagrams , including the one from the Z penguin, are 
negligible. In this case we have 

4 

= 1.02 ■ 10"i\ (81) 

which is independent of b. 

The range of values for BR{Kl fie) agrees with the actual experimental upper bound 
of 3 ■ 10~^^ at 90%CL [|1^. This range is however quite wide, spanning eight orders of 
magnitude, and is due to the fact that the branching ratio has a strong dependence on the 
masses of the vector-like fermions F''"? and Aii"^. Notice from Table 1 that F''"? are the 
vector-like leptons and quarks which are SU{2)l doublets, while jWf denotes an SU{2)l 
singlet charged lepton and quark respectively. It is interesting to note from Fig. 10 that, 
in the case where all vector-like fermions are in the interesting mass range (a few hundreds 
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— lOOTeF 



, 327r4M|!os2 ^Q^fC,{xF,)Co{xF,] 



BR= 11.24- 10"^2 



gp lOOTeV 
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M, 



G 



of GeV's) where they can be produced by QCD or Electroweak processes, BR{Kl /ie) 
is hopelessly small to have a chance to be observed (the lower fiat part of Fig. 10). At the 
other extreme (the upper flat part of Fig. 10), BR{Kl ~^ /ue) ~ 10~^^ while the mass of 
TWi"^ is unreachable (of 0(few hundreds TeV's) assuming the mass of F''"? is of order a few 
hundreds of GeV's) by any conceivable earthbound machine. 

IX. EPILOGUE 

We have presented a link between the quark sector and the lepton sector based on a model 
of neutrino masses 0. This link manifests itself in a common (small) parameter, 6, which 
appears in both sectors. In the neutrino sector, h is crucial in determining the magnitude 
of the mass splitting Am?, and, in an indirect manner, the masses themselves. In fact, the 
model presented in Ref. deals with the case of three light, degenerate neutrinos, and the 
lifting of this degeneracy is proportional to h. In the quark sector, as we have seen above, this 
same parameter appears in various FCNC penguin operators which are seen to be linear in h 
at the lowest order. It is well-known that these penguin operators are important in various 
aspects of Kaon physics: e'/e, rare K decays, etc... Strong constraints in this sector would 
also constrain the neutrino sector as well. It turns out, as we have shown above, that there 
is no contribution to e'/e from our model, to lowest order in h. However, at the same order 
in 6, our model makes a contribution to the rare decay process — > Taking into 

account the upper bound on the short distance contribution to that process, we found a 
strong constraint on h which is 6 < 2.5 x 10^^. 

We have also calculated the BR{Kl /le), finding that it strongly depends on the 
masses of the vector-like fermions F''''^ and A^i^ and on b for a certain range of those masses. 
We found that BR{Kl fie) goes from 6.65 ■ 10^^^ which makes this decay practically 
unobservable, to 1.02- 10"^^, with the choice of Mq = 200TeV^ and xp = 10^^, when xm/xf 
goes from 1 to 10^^. For b we have used the upper bound derived from fi^fx^. As 

one can see from Fig. 10, the first case (6.65 ■ 10^^^) corresponds to the interesting situation 
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where it might be possible to produce and observe these vector-hke particles since their 
masses could lie in the few- hundred GeV region |T^, despite the fact that one will not be 
able to observe Kl — > yue. The second case (1. 02-10"^^) corresponds to a possible observation 
for Kl — > /ie, while forsaking that of the singlet quark and charged lepton jWi^. And finally, 
there are these in-between cases as can be seen from Fig. 10. 

Because of the upper limit on 6, and obviously because of the choice of heavy family 
gauge bosons, Mq of C(lOOTeV^), it turns out that the physics of the kaon sector by itself, 
in our model, is not too different from the SM. Nevertheless, we have seen that there is 
still some margin for possible contributions of new physics to Kl — >■ n~^fi~, the bound on 
K^ — s> TT^uu. As for the decay Kl fJ^e, which is forbidden in the SM, we have seen 
that our model can make a non-negligible contribution {BR ~ 10^^^) which is practically 
independent of b. In the region where it depends strongly on b (the lower flat region of Fig. 
10), the branching ratio is negligible, practically similar to the SM with a Dirac neutrino. 
In that region, as we have stressed above, the new physics signal would be the production 
and observation of the vector-like fermions. 

The bound on b could have interesting implications on the neutrino sector itself, as we 
have mentioned above. It was shown in Ref. how the parameter b affects the mass 
splitting of three formely degenerate neutrinos. In particular, it was shown how Am^ is 
sensitive to 6. However, it was also shown how b indirectly affects the overall magnitude 
of the masses. Unfortunately, at the present time, one is quite far experimentally from a 
direct determination of the masses themselves. Needless to say, future experiments are of 
paramount importance to this crucial question. 
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X. APPENDIX A 



In this appendix, we will summarize the results of Ref. 0. The Yukawa part of the 
Lagrangian involving leptons can be written as 

^lepton = 9Ell(j)eaR + GJin^Fn + GM^FL(j)MiR + Gm2Fl4>M2R + GiMiL^aeR + 

G'3-M2lP° Cfl + MfFlFr + MmMilMir + Mm,M2lM2r + h.c. . (82) 

which, after integrating out the vector-like fermions F, Aii and Ai2, brings to the effective 
Lagrangian 

GMll{nJp^)rfpR + h.c., (83) 

where 

_ G1GM1G2 ^ _ G1GM2G3 ,r,A\ 

- M,Mm, ' - M,Mm, ■ ^^^^ 
The main assumption in building the above Lagrangian is the conservation of lepton 
number L, thereby forbidding the presence of Majorana mass terms. The way < Vt >= 

(i\ 

(0, 0, 0, V) and < p >= (0, 0, 0, V ® Si), with si = 



have been chosen, makes the 



vV 



neutrino of the 4*^ family massive at tree level, while the others three neutrino remain 
massless. These three neutrinos would get a mass dynamically via loop diagrams. For the 
neutrino of the 4*^^ family we have 

^n = G,Gm,G,j^^^^, (85) 

which, for 



VV'/MfMm2-0{1), (86) 

can be expected to be even as heavy as 175GeV, satisfying the bound of Mz/2 from LEP. 
For the three light neutrinos, one has 
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m, = mN^^, ^^^I{G,P), (87) 

where A/(G', P) is given by 

1 M^[M|(M^ln(£|)-M|.ln(£|)) + M^M2ln(g|)] 

(88) 

The main result that was obtained in Ref. ||^ is that one can obtain for rriy a value of the 
0{eV), or equivalently the ratio R = rriu/mN ~ O{10~^^) as long as the ratios of masses of 
the particles propagating in the loop satisfy certain relations. In Ref. |^ it has been shown 
that, taking the masses in units of Mp, where Mp could be of 0{> 200GeV), one obtains 
R < 10-" for Mg/Mm2 ^ 10~^ when Mg < 10^ or for Mg/Mm^ ~ 10"^ - lO^^ when 
Mg > 10^, with Mp ~ 1 — 10^. After lifting the degeneracy of the three light neutrinos, by 
breaking the remaining family symmetry 5*0(3), one obtains the following mass eigenvalues 
for the three light neutrinos 

mi = m/-^^{AI{G, P) - bAI{G, P,b)} , (89) 
m2 = mM'-^^AI{G,P), (90) 



ma = m/-^^{AI{G, P) + 6A/(G, P, -6)} , (91) 
oZ 



where AI{G, P, ±6) is given by 

AI{G, P, ±6) = I{Mg, ±b) - I{Mp, ±b) , (92) 

with 

Mp[-MUl±b + H^)) + MUl±b)] 

(93) 

The mass splittings, neglecting terms of order b^ are now 
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ml-ml = m2mN2b^-^^^AIiG, P, -h) , (94) 

ml-ml = ,n2m;v26^^^A/(G, P, h) . (95) 

The above mass splittings are almost degenerate and in Ref. 0] a possible solution to lift this 
degeneracy was presented. What has been showed is that if one introduces mixing terms in 
the neutrino mass matrix between the three light families and the 4*'' family, one can obtain 
reasonable values for the mass splittings. Two numerical examples have been presented. 
The mass splittings for the particular choice of the parameter h = 0.000095 that satisfies 
the upper constraint for b derived in this paper and presented in Eq. ^ are given in Eq. ^ 
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XI. APPENDIX B 

In the following we will show in a simple example how the dependence on the parameter 
b, introduced in Eq. |16|, for the amplitude of the FCNC processes in our model appears. 
This example will be introduced ill cLll cLIIclIcL gous fashion to a FCNC process in the SM, 
which we know to be governed by the GIM mechanism. The process that we are considering 
is t — + cW^W~ . Now the amplitude of this process is proportional to the quantity 

j=d,s,b '"'3 

where the sum appears because one has to consider all possible virtual down-quarks, and p 
is the momentum of the virtual quarks. Notice that Eq. ^ represents the approximation 
p ^ nij which is a valid one as we shall see below. Vij is the CKM matrix and the unitarity 
of this matrix implies in particular the relation 

E ^t;-K. = 0. (97) 

j=d,s,b 



Using the above relation we can write Eq. 96 as 
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v*v 



+ 



2 2 

rrij — 



EV*V 



(98) 



Using now the approximation 



(99) 



and doing the sum we can rewrite Eq. as 



mt — mz 



From kinematic considerations falls in the range 



2\ • 



(100) 



{Mw + m^f < p2 < (rn, - Mwf , 



(101) 



which means that p^ is of order of M^. One obtains the final expression 



e 



2 "^6 1 



(102) 



which shows the characteristic term of the GIM mechanism ml/M^ for FCNC processes 
involving external t. 

Now, looking at our model, and assume that the process t cF A4i can kinematically 
occur, (although it cannot be in reality because Mp > 200GeV and M_m^ > Mp). The 
amplitude will be proportional to 

1 



T 



(103) 



where Aq^ is a unitary matrix given by (see Ref. 



Ac 



( I ^ 1 \ 

2 ^2 2 

1 1_ 1 

2 V2 2 

1 n 1 



(104) 



chosen different from that one given in Eq. |^ because otherwise we will have no transition 
from top to charm, and the sum comes out from the fact that we have to consider as virtual 
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states any gauge boson of 5*0(3). The unitarity of the matrix imphes the particular 
relation 

E^nj3<2, =0, (105) 

j 

that we can use to rewrite the amplitude as 

- + ^ A.^^.,Al,,^^^ . (106) 

in the same way as we did for the SM. Now for the particular matrix Aq that we have 
chosen, we obtain 

1 



AnMi23 = , (107) 

which gives 



1 - 

2V2{k'-Ml){k^-Ml)- 

Substituting for Mq^ and Mq^ the expressions given in Eq. |1^ we obtain 



1 bM^ 



v/2 (F _ M2 (1 + 6))(P - M2 (1 - b)) ■ 
From kinematically considerations we have the following constraints for k'^ 



(109) 



(m, + MmJ' <k^< {mt - Mpf , (110) 



;iii) 



which implies /c^ <C M^, with M^^, Mp <C Mq. The final expression becomes 

1 b 
~^k^-Ml ' 

which shows the linear dependence on the parameter b. One notices that the linear depen- 
dance on the parameter b comes out in the same way as the linear dependence on the quantity 
m^/M^ in the GIM mechanism of the SM. It is important to notice here that, in the SM, 
the GIM suppression comes from two factors: the appropriate product of CKM elements 
and the ratio m^/m^, where is a quark mass (up or down). (Even if m^/m^ ~ C(l) 
such as for the top quark, the CKM suppression factor can be very small.) In our case, the 
amplitude is always suppressed by the parameter b. 
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FIGURES 

FIG. 1. Feynman graphs for the effective verteces of the Z and 7 for AS = 1 processes s ^ d. 
The self energy diagrams are omitted. 

FIG. 2. Feynman graphs for the effective vertex of the g for AS = 1 processes s ^ d. The self 
energy diagrams are omitted. 

FIG. 3. The function Co{xF,b), appearing in the amplitude for the Z effective vertex is plotted 
in linear and logarithmic scale for b = 0.001 and b = 0.0001. The function Co{xf), which is the 
first non vanishing coefficient in the Taylor expansion in 6 of Co (a::^! b) is also plotted. 

FIG. 4. The logarithm of the function Do{xF,b), appearing in the amplitude for the 7 effective 
vertex is plotted versus the logarithm of a;^ for 6 = 0.001 and b = 0.0001. The function Dq{xf), 
which is the first non vanishing coefficient in the Taylor expansion in b of Do{xf, b) is also plotted. 

FIG. 5. Box diagrams for AS = 1 process sd dd. The corresponding box diagrams with the 
pseudo NG bosons Rep'^ instead of the NG bosons are omitted. 

FIG. 6. The function Bq{xf, b), appearing in the amplitude for the box diagrams in the process 
sd dd is plotted in linear and logarithmic scale for b = 0.001 and b = 0.0001. The function 
Bo{xf), which is the first non vanishing coefficient in the Taylor expansion in b of Bo{xf, b) is also 
plotted. 

FIG. 7. The function Bmixo{xF, xm, b), appearing in the amplitude for the box diagrams in the 
process sd dd is plotted in linear and logarithmic scale for 5 = 0.001 and b = 0.0001 versus xm, 
with the condition xp = 10~^a;M- 

FIG. 8. Feynman graphs for the AS = 1 process sd — /xe. In the graph (e) the vertices are 
effective. The corresponding diagrams with the pseudo NG bosons Rep'j^ instead of the NG bosons 
J^'j are omitted. 
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FIG. 9. The logarithm of the function {CLg — Cjiq)'^{Cli + Cjy^), which appears in Eq. 76 for the 
BR{Kl lie) is plotted versus Log{xM/xF), for 6 = 2.5 • 10"^, Mq = 200TeV and xp = 10"^. 

FIG. 10. The logarithm of BR{Kl fJ-e) is plotted versus Log{xM/xF), for h = 2.5 • 10~^, 
Mg = 2mTeV andxF = 10"'^. The lowest flat zone corresponds to BR{Kl fie) = 6.65 • lO"^^, 
while the highest flat zone corresponds to BR{Kl — > fie) = 1.02 • 10"^^. 
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TABLES 

TABLE I. Particle content and quantum numbers of 

SU{3)c ® SU{2)l ® U{1)y ® SO{Nf ) ® SU{2)^^ 









/r — fl 2 —1 /2 Nt 1) 




r>, iQ ^ O /Q AT 1\ 




— ('^ 1 —1 AT^X Tl 
"il — I, -J) J-, -L/o, iv/, Ij 






Right-handed z^'s 


Option 1: r]R= {l,l,0,I\f,2) 




Ontinn 9- no — MID 9V 




„/ _ fi 1 n 1 0) 


Vector-hke Fermions 


-^L,i? = (li 2, -1/2, 1,1) 
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